Equation of state of a multicomponent G?-dimensional hard-sphere fluid 
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A simple recipe to derive the compressibility factor of a multicomponent mixture of d- dimensional 
additive hard spheres in terms of that of the one-component system is proposed. The recipe is based 
(i) on an exact condition that has to be satisfied in the special limit where one of the components 
corresponds to point particles; and (ii) on the form of the radial distribution functions at contact 
as obtained from the Percus-Yevick equation in the three-dimensional system. The proposal is 
examined for hard discs and hard spheres by comparison with well-known equations of state for 
these systems and with simulation data. In the special case of d = 3, our extension to mixtures of 
the Carnahan-Starling equation of state yields a better agreement with simulation than the already 
accurate Boublfk-Mansoori-Carnahan-Starling-Leland equation of state. 



Due to their importance in liquid state theory, for years 
researchers have proposed empirical or semi-empirical 
(analytical) equations of state of various degrees of com- 
plexity for one-component hard-sphere fluids. Notable 
among these, is the celebrated Carnahan-Starling (CS) 
equation of state |Q , which is not only rather simple but 
also accurate in comparison with computer simulation 
data. In the case of hard-sphere mixtures, the proposals, 
also empirical or semi-empirical in nature, are much more 
limited, with the Boublfk-Mansoori-Carnahan-Starling- 
Leland (BMCSL) equation Q standing out as the usual 
favorite. The situation for one-component hard-disc flu- 
ids is rather similar. Here, no analog of the CS equation 
using the |(v) + f(c) recipe has been derived, due to the 
absence of an analytical solution of the Percus-Yevick 
(PY) equation in this instance. Nevertheless, accurate 
and simple equations of state have been proposed, such 
as the popular Henderson equation and the recent 
one by the present authors Hard-disc mixtures, on 
the other hand, have received much less attention and the 
proposed equations of state for these systems are rather 
scarce. Given this scenario, the major aim of this paper 
is to show that, on the basis of a simple recipe, accurate 
equations of state of a multicomponent d-dimensional ad- 
ditive hard-sphere mixture may be derived, requiring the 
equation of state of the one-component system as the 
only input. The recipe makes use of a consistency con- 
dition that arises in the case that one of the components 
in the mixture has a vanishing size, as well as from some 
insight gained from the form of the radial distribution 
functions at contact given by the solution of the Percus- 
Yevick equation for a hard-sphere fluid in three dimen- 
sions ||. 

Let us consider an Y-component system of hard 
spheres in d dimensions. The total number density is 
p, the set of molar fractions is {x\, . . . , xn}, and the set 
of diameters is {a\, . . . , <jn}- The volume packing frac- 
tion is 77 = v d p{a d ), where v d = (tt / %) d l 2 /r(l + d/2) is 
the volume of a d-dimcnsional sphere of unit diameter 
and (a n ) = JY Xicr™. In the case of a polydisperse mix- 
ture (Y — > oo) characterized by a size distribution /(<r), 



(a n ) = J dacr n f(a). Our goal is to propose a simple 
equation of state (EOS) for the mixture, Z^ N >{rj), con- 
sistent with a given EOS for a one-component system, 
Z^irf), where Z = p/pkuT is the usual compressibil- 
ity factor. A consistency condition appears when one of 
the species, say the Y-th, has a vanishing diameter, i.e. 
<7jv — * 0. In that case, 
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At a more fundamental level, we will consider the contact 
values of the radial distribution functions, g\f\<Jij), the 
knowledge of which implies that of the EOS through the 
relation 
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Taking as a guide the form of obtained through 
the solution of the PY equation for hard spheres (d = 3) 
H, we propose to approximate by a linear interpo- 
lation between and (1 — r/)^ 1 , namely 
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When the above ansatz is inserted into equation 
gets 
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This form of the EOS complies with the requirement ([l]). 
Note that equation (0) expresses Z^ N >{rj) — 1 as a lin- 
ear combination of Z^'(rf) — 1 and (1 — 77) _1 — 1, but 
the dependence of the coefficients on the size distribu- 
tion is much more involved than in equation ([}]). The 
key outcome of this paper is the EOS given by equation 
in which the compressibility factor of the mixture 
is obtained from that of the one-component system for 
arbitrary values of the dimensionality d and the num- 
ber of components N. It is worth noticing that in the 
one-dimensional case, equation (0) yields the exact re- 
sult (77) = Z« (77). 

As a straightforward application of equa- 



tion 



the virial 



one 
,i-d, 



can 
(i) 



2 i -< 1 Ao&: 

coefficients 

71 = 2 



easily 
1- A 



get 



B 



where 



defined by 
(i) 



Bn ' are 

1 + ZrT=2 B « V) p™" 1 - and where U 
are the reduced virial coefficients of the one-component 

system, i.e. ZW(rj) = 1 + En= 2 &£V _1 . 

We will now focus on the case of hard discs (d = 2). 
Equation (Q) then becomes 
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The relationship between Z^ N '(rf) and Z^ 1 ^) as given 
by equation (^|) rests on a different rationale from that 
pertaining to another simple proposal, namely the Con- 
formal Solution Theory (CST) j|, 0. In this latter the- 
ory, the EOS reads Z^\,(r/) = Z^(rj c g) with r) e a = 
i (l + (cr) 2 /(er 2 )) 77, but we note that this equation does 
not comply with the general requirement (jl]) . If the one- 
component system is assumed to be described by the 
Scaled Particle Theory (SPT) our extension to mix- 
tures takes on a particularly simple form: 
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This is precisely the true SPT EOS for mixtures 0, 
which is indeed rewarding. The EOS Z^ (77) = 

[l — 277 + (2?7o — 1)?7 2 /?7q] , where 770 = V3ir/6 is the 
crystalline close-packing fraction, has been recently pro- 
posed by us 0] to describe a one-component system. 
When this EOS, hereafter referred to as the SHY EOS 
following the nomenclature introduced in reference ||, 
is substituted into equation we obtain the following 
extension: 
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The well-known Henderson (H) EOS |3| can also be ex- 
tended: 

w _ l-(l-(a) 2 /( a 2 )) ,7 + (6« - 3) ((a) 2 /^ 2 )) t? 2 
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-\/3- This equation is quite similar 
to the one proposed by Barrat et al. 0, the only dif- 
ference being that the coefficient of ?7 2 in the numerator 
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where 63 is the exact (re- 
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duced) third virial coefficient. Nevertheless, since 63 is 
well approximated by 1 + (b^ - l)(cr) 2 /(cr 2 ), both EOS 
are practically indistinguishable. Although we could con- 
sider the extensions of other EOS originally proposed for 
a one-component system of hard discs (for a list of many 
such EOS we refer the reader to references [|), for 
the sake of simplicity we will restrict our analysis to the 
SPT, eSHY and eH EOS. 

Let us now consider the virial coefficients Sn for 
hard discs. It follows that in this 
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tion yields the exact second virial coefficient [Q, the 
higher coefficients being approximate. In the par- 
ticular case of a binary mixture, the composition- 
independent coefficients Bn},n 2 ar e defined through 

According 
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where 71 = 711+712 and a — o~ilo~\. This form has the 
same structure as the interpolation formula suggested by 
Wheatley jn^. In fact, he proposes an EOS (henceforth 
labelled as W) of the form 
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where the coefficients c„ are chosen so as to reproduce 
the first eight virial coefficients given by the interpolation 
formula |Tc[ . 

Now, let us consider the case d = 3. Equation (^) then 
yields 
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Using the CS EOS §, Z«(r?) = (1 +77 + 7 7 2 — r/ 3 )/ ( 1 — r/) 3 , 
the result may be expressed as 
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where the BMCSL EOS § is 
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As another, example, let us consider the Carnahan- 



Starling-Kolafa (CSK) EOS Z^{rj) 
2rf{\ + 77)73] /(l - ?7) 3 . Its extension is 
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which docs not coincide with Boubh'k's extension to mix- 



tures of the CSK EOS Jl2 
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Recently, Henderson and Chan (HC) have proposed a 
modification of the BMCSL EOS [TJ, [TJ for the partic- 
ular case of a binary mixture in which the concentration 
of the large spheres is exceedingly small, starting from an 
asymmetric prescription for the radial distribution func- 
tions at contact. The resulting EOS, with a\ > 02j is 
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We shall now perform a comparison with the (very 
few) available computer simulation data. We begin with 
hard-disc mixtures. In figure ^ we display the packing- 
fraction dependence of the compressibility factor Z for 
the SPT, W, eH and eSHY EOS, together with the sim- 
ulation results of Barrat et al. |7| , for the binary mixture 
defined by x% — 0.351 and o%J(T\ = 0-8- In this case, 
the performance of the eSHY EOS is outstanding and 
clearly superior to all the other choices. To complete the 
picture, in figure |2] we present the results for the ratio of 
the fifth virial coefficient to the fourth power of the (ex- 
act) second virial coefficient as a function of the larger 
disc concentration and for two size ratios. Here, the best 
agreement with the numerical data of Wheatley |TjJ is 
obtained with the eH EOS, which is not very surpris- 
ing since in the one-component case (x± = 1) it gives a 
very good estimate of this ratio. Nevertheless, the over- 
all trends including the position of the maximum are still 
captured in all approximations. 

As far as hard-sphere mixtures are concerned, the 
following comments can be made. To our knowledge, 
only simulation results for binary mixtures have been 
reported. The most recent data Il6|] indicate that the 
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FIG. 1: Compressibility factor as a function of the packing 
fraction for a binary mixture of hard discs defined by x\ = 
0.351 and cr 2 /<Ji = 0.8. 
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FIG. 2: Fifth virial coefficient as a function of 11 for two 
binary mixtures of hard discs defined by o^/oi =0.1 and 
a 2 /ai = 0.2. 



BMCSL EOS underestimates the pressure as obtained 
through simulation. In fact, the BCSK EOS is geared to 
correct this deficiency, at least for 77 < 0.5, in a similar 
fashion as the CSK EOS corrects the CS EOS. As a gen- 
eral trend, for 77 < 0.5, our extended EOS, namely the 
eCS and the eCSK, also go in the correct direction. More- 
over, in this density range, Z B mcsl < Z cC s < Z oC sk and 
Zbmcsl < Zbcsk < Z cC sk- Although limited in scope, 
the results shown in figure ^ illustrate these features. 
Here we have considered an equimolar binary mixture 
with size ratio 02 /ai = 0.6. As the differences between 
the values predicted by the various EOS for the com- 
pressibility factor Z are very small, we have chosen to 
present the results, including the simulation data of Yau 
et al. E3 (open circles) and Barosova et al. [Il6| (filled 
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FIG. 3: Compressibility factor as a function of the packing 
fraction, relative to the BMCSL value, for an equimolar bi- 
nary mixture of hard spheres with 02 /or = 0.6. 
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FIG. 4: Fifth virial coefficient as a function of xx for two 
binary mixtures of hard spheres defined by as /ax = 0.05 and 
ail ax = 0.2. 
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circles), in terms of the packing fraction dependence of 
Z — Zbmcsl- Despite the scatter of the simulation re- 
sults of Yau et al. [ p^[ , it is apparent that, depending 
on the range, both the eCSK and eCS EOS seem to do 
a better job than either the BMCSL, BCSK or HC EOS 
(although in all fairness we should add that the equimolar 
condition is beyond the scope for which the latter EOS 
was originally devised). In fact, if one considers a higher 
density point computed by Yau et al. |l4j (rj — 0.59, 
■^simui — -Zbmcsl = 0.26) that is off the scale, it appears 
that the overall trend is better captured by the eCS EOS, 
although for rj < 0.5 the eCSK should probably be the 
preferred EOS. The results for the composition depen- 
dence of the fifth virial coefficient for a binary mixture 
and two values of oi / <j\ displayed in figure ^ also indi- 
cate that all the approximations lead to very good val- 
ues as compared to the recent numerical data of Enciso 
et al. with a slight superiority of the BCSK EOS 

for the region around the maximum. The HC and the 
eCSK results have not been included, since they are al- 
most identical to the ones of the BMCSL and the eCS, 
respectively. 

In conclusion, it is fair to state that we have introduced 
a very simple and general recipe that allows one to get 
a reasonably accurate approximation to the EOS of a 
multicomponent mixture of d-dimensional hard-spheres 
from any reasonable EOS of the one-component system. 
It also seems that, as exemplified in the case of binary 
three-dimensional hard-sphere mixtures, the more accu- 
rate the EOS of the one-component system, the better 
results the approximation yields. 
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